Response of a d x 2_ y 2 Superconductor to a Zeeman Magnetic Field 



Kun Yang 1 and S. L. Sondhi 1 * 2 

1 Department of Physics, Princeton University, Princeton, N J 08544 
2 Institute for Theoretical Physics, University of California at Santa Barbara, Santa Barbara CA 93106-4030 

(February 1, 2008) 

We study the response of a two dimensional d x 2_ y 2 superconductor to a magnetic field that couples 
only to the spins of the electrons. In contrast to the s-wave case, the d x 2_ y 2 state is modified even at 
small magnetic fields, with the gap nodes widening into normal, spin polarized, pockets. We discuss 
the promising prospects for observing this in the cuprate superconductors in fields parallel to the 
Cu-0 planes. We also discuss the phase diagram, inclusive of a finite momentum pairing state with 
a novel linkage between the momentum of the pairs and the nodes of the relative wave function. 



Following the original work of Clogston and Chan- 
drasekhar fw, the modification of superconductivity by 
the Zeeman coupling between the spins of the electrons 
and an applied magnetic field has attracted intermittent 
attention ||] . Much of this has centered on the bound on 
the upper critical field H C 2 provided by consideration of 
the Zeeman interaction alone ( "Pauli limit" ) and on the 
nature of the phase boundary when this effect dominates. 
A more unusual aspect of this physics was uncovered by 
Fulde and Ferrell || and by Larkin and Ovchinnikov j|] 
in the possibility of a finite momentum pairing state at 
large magnetic fields, where the reference Fermi surface 
is spin-split. Experimentally, the classic predictions on 
the nature of the s-wave phase boundary have been borne 
out by work on thin Al films Q , while recent work on un- 
conventional superconductors has exhibited Pauli limited 
critical fields || as well as the first evidence for a real- 
ization of the finite momentum pairing state j7j. Very 
recently, the Zeeman suppression has also been discussed 
in experiments on mesoscopic samples ||. 

In this paper, motivated by the recent experimental 
identification of the pairing state in several of the cuprate 
superconductors, we discuss the Zeeman response of an 
ideal two-dimensional d x 2_ y 2 superconductor. We be- 
lieve that this is a useful exercise on three grounds. First, 
the cuprates are strongly two-dimensional systems and 
hence admit a geometry for measurements in a mag- 
netic field, with the field parallel to the Cu-0 planes, 
where the Zeeman response should dominate the orbital 
response at low temperatures. Second, we estimate that 
spin-orbit scattering, which attenuates the Zeeman re- 
sponse, is small enough in the cleanest samples to allow 
its neglect above fields as small as a few tesla. Third 
and most interestingly, the existence of nodes in the gap 
function imply that (in contrast to the s-wave case) the 
superconducting state responds non-trivially at arbitrar- 
ily small values of the magnetic field. As is intuitively 
plausible, the response is paramagnetic with the destruc- 
tion of superconductivity over parts of the Fermi surface 
where the Zeeman energy \iB exceeds the local gap A(k), 
and a spin polarization of the resulting normal electrons. 



This leads to observable effects in all quantities that are 
sensitive to the density of states for quasi-particles. 

In the following we will explicitly illustrate this by cal- 
culations on a weak coupling BCS model and present es- 
timates that indicate that it is readily observable in the 
cuprates. For completeness, we also discuss the mean- 
field phase diagram of our model in the field-temperature 
plane, where we note the existence of finite momentum 
pairing with a novel linkage of the wave vector to the 
nodal structure of the gap function. 
Choice of Hamiltonian: We study a two-dimensional 
(2D) electron system described by the Hamiltonian 

H= efc4 CT c k^ + X) Fk ' k ' c k,T c -k,l c -kM c k',T' (!) 

k,o-=T4 k,k' 

where is the rotationally invariant kinetic energy mea- 
sured from the Fermi energy e^, and for |efej, \e' k \ < e c <C 
6f, the pairing potential V takes the form |^] 

Vt, k ' = -2V cos(20 fe ) cos(20fc/) + /i*, (2) 

where Vq > 0, /i* is the renormalized s-wave repulsion 
and 9k is the azimuthal angle of k; V = otherwise. At 
low temperatures the system is a superconductor with a 
gap function of the d x 2_ y 2 form: A(9) = Ao(T) cos(2#), 
where Ao satisfies 

N(0)Vq [ 2 * r° 2 cos 2 (29) , ( m9) 

and E(£,9) = ^ 2 + A 2 cos 2 (2(9), N(0) is the single 
particle density of states for each spin species at the 
Fermi level. In weak coupling N(0)Vq <C I, assumed 
throughout this paper, this leads to a maximum gap 
A 00 = A (T = 0) = 2.43e c e- 1 / JV (°) v ' , and a critical tem- 
perature T c = 0.467Aoo- The quasi-particle spectrum is 
governed by the mean-field Hamiltonian: 

H MF = £k(4"k + 0tpk), (4) 

k 

where = ^/e 2 (k) + Aq cos 2 (2#fc), and a k (/?£) are cre- 
ation operators of up (down) spin quasi-particles. 
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We should note that the d x 2_ y 2 state that we are in- 
terested in, arises in (nearly) tetragonal lattice systems. 
Our choice of model here is intended to mimic this lat- 
tice physics with minimal and computationally favorable 
ingredients: we keep a rotationally invariant Fermi sur- 
face but introduce a pairing potential (^) that has the 
reduced symmetry of the lattice. 

A uniform magnetic field (B) applied parallel to the 
2D plane does not couple to the orbital motion of the 
electrons in the plane. It does, however, lift the spin 
degeneracy, and introduce the Zeeman term, Hz = 
— fiB X)k( c kT Ck T — c iki Ck l)' t ne Hamiltonian, where 
(i = gy, B /2 is the magnetic moment of the electron, and 
t and I refer to spin direction along and opposite to the 
field direction respectively. This modifies both the gap 
equation and the quasi-particle Hamiltonian to 

i = »m f : de r v os2(2 * } 
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X [tanh( 2k B T )+tanh( 2k B T )] 
H MF = - MS)4«k + (Ek + nB)0l0k)] ■ (5) 

k 

Qualitatively, the Zeeman field lowers/increases the en- 
ergy of the spin up/down quasiparticle states which 
in turn changes their occupation and affects the self- 
consistency condition for the gap function (which is now 
distinct from the true quasiparticle gap). 
Weak-field Response: We begin by considering weak 
magnetic fields and low temperatures: [iB, k B T <C Aoo- 
An s-wave state is essentially unaffected in this limit. 
This is because the occupation numbers for quasiparticle 
states remain exponentially small at low T and B, due 
to the finite gap, even though the field shifts the quasi- 
particle bands linearly (at T — the gap function and 
the ground state are completely unaffected). 

The situation, however, is qualitatively different in the 
case of the d x ?_ y i state studied here; in our case the gap 
vanishes at four nodal points on the Fermi surface; suf- 
ficiently close to these points there are always some spin 
up quasiparticle states whose energies become negative 
for arbitrarily small values of B. These states, which live 
in elliptical pockets, < fj,B, near the Fermi surface 
(Fig. 0), develop a thermal occupation that is 0(1) at 
any temperature. Translating back into electron opera- 
tors, one sees that these pockets are in fact regions of the 
spin polarized normal state (fully polarized at T = 0) — 
their inner and outer arcs are pieces of the spin split 
Fermi surface which come together when the angle de- 
pendent gap function exceeds the Zeeman energy. The 
loss of pairing in them and the area within k B T, leads to 
an overall reduction of the gap function: 

A (T, B) = A 00 [l - (^) 3 F A (-^-) + O(A " 4 )] 



Fa(x) 



t 2 [l - (tanh 



t + x 



tanh 



t — x 



)/2}dt . (6) 



Noting that -Fa(O) « 3.61 recovers the zero field answer, 
while for x 3> 1 Fa(x) ~ | a; | 3 /3 whence, 



A (T = 0,B) = A 00 [l - l(^) 3 + 0(S 4 )] 



r An 



(7) 



Consequently, the reduction of the gap function at low 
fields and temperatures is quite modest, and the most 
important effect of the proliferation of quasi-particles in 
the ground state is the finite density of states (DOS) at 
the Fermi level, which qualitatively alters the low energy 
physics of the system [ll]]. We now turn to the conse- 
quences for some physical quantities in this regime, where 
we may neglect the reduction of An. at leading order. 
Specific Heat: This takes the scaling form 



C{T,B)=2k B N{Q) { ^^Fc{^=) 



Fc(x) = E / *(* 



A o k B T' 

2 t+ax /f t+ax , -\\2, 



(8) 



+ axYe t+rTX /(e t+ax + iydt 



For fields in excess of the temperature 



2g-a; whence 



x > 1, F c (x) ~ 
C = 2fk B TN(0)fxB/A 00 = C N[ LBj A 00 , 
where Cat is the normal state specific heat; the linear T 
dependence at low T is a consequence of the finite DOS. 
For iiB <§; k B T, we recover the expected T 2 dependence 
upon using F c (0) = 9C(3) w 10.8. 

Thermal Conductivity: At low temperatures where im- 
purity scattering may lead to a constant quasiparticle 
scattering rate [12 ], the thermal conductivity n e should 
be proportional to C . Therefore n e should increase lin- 
early with T, for k B T <C ^B, while it increases with T 
quadratically in the absence of the field, as observed ex- 
perimentally [[u| . This is to be contrasted with the recent 
experimental finding jl3| that a magnetic field applied 
perpendicular to the Cu-0 plane of the cuprate super- 
conductor suppresses the electronic thermal conductivity, 
presumably due to orbital effects. 
Magnetization: M also takes a scaling form: 

M(T,B) = 2vN(0) { -^p^F M A 
A o k B T 



F M (x) 



1 



1 



1 e t+x + 1 



)dt. 



(9) 



^k n T' 



The limits F M (x > 1) - x 2 /2, and F M (x < 1) ~ 
(2 log 2)a;, imply M cx B 2 when fiB > k B T, and M cx B 
when fiB <C k B T. 

Tunneling Conductance: The T = tunneling conduc- 
tance of a superconductor-insulator-metal junction at 
varying bias G(V) is, in principle, the most direct mea- 
surement of the single particle DOS of the superconduc- 
tor. G goes to zero linearly with V for a d x 2_ y 2 supercon- 
ductor in zero field. For fiB ^> k B T, the finite DOS leads 
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to a finite conductance: G(V — 0) = G n ^B/ Aoo, where 
G n is the tunneling conductance when the superconduc- 
tor is in its normal state. For eV < fiB, the tunneling 
current is spin-polarized |l4[ . 

Phase Diagram: We now turn to the effects of strong 
magnetic fields. An s-wave superconductor undergoes 
a first order phase transition to the normal state when 
fiB = A/ \/2 (A is the s-wave gap) at T = Q, ignoring 
the finite momentum pairing instability (see below). This 
follows upon noting that the singlet s-wave state is insen- 
sitive to the Zeeman field, while the normal state lowers 
its energy in proportion to its Pauli susceptibility. The 
temperature-magnetic field phase diagram exhibits a 
tricritical point where the first order line terminates, and 
the field tuned transition becomes continuous. 



The cL 



superconductor does respond to the Zee- 



man field but far more weakly than the normal state. 
This leads to a phase diagram (Fig. 2) of the same topol- 
ogy as in the s-wave case. We find that at T = there is a 
first order transition to the normal state at \xB = 0.56Aoo 
(close to the value 0.5Aoo obtained without account- 
ing for the paramagnetism of the d x i_ y 2 state); at the 
transition the gap function has amplitude 0.92 Aoo- For 
T < 0.56T C , the transition remains first order, and the 
normal (superconducting) state becomes local minimum 
of free energy at lower (higher) temperature, as repre- 
sented by the dotted lines in Fig. 2; above T — 0.56T C 
the transition becomes continuous. 

Thus far we have only considered zero momentum pair- 
ing. At high fields one may suspect that it might be 
more favorable to try to pair across the spin-split Fermi 
surface of the partially polarized normal state; indeed 
it is known in the s-wave case j^,|],|l0| that this hap- 
pens at low temperatures. The pairing is then at a fi- 
nite center of mass momentum for the Cooper pairs and 
the resulting transition from the high field normal state 
is continuous. At lower fields this Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state gives way to the zero momen- 
tum pairing state by a first order transition. 

In order to consider this possibility for the d x 2_ y 2 
problem, we extend the pairing potential in ([!]) to al- 
low pairing between electrons with total momentum q, 
for q <C e c /vp, where vf is the Fermi velocity: 

V= X] ^ k , k <4 T cl k+qa c_ k < +qa c k , !T . (10) 

k,k\q 

As we will see later, the pairing momentum q is at most of 
order A o/vf, which is much less than € c /vf- Therefore 
it is a reasonable approximation to neglect the depen- 
dence of the pairing matrix element Vic,k' on q. The gap 
still takes the d x 2_ y 2 form and obeys 

i=^mr d e fv os2( ^ 



2?r 



E(JZ, 



, , ,J5(f,0) +z(6). , ,E(£,6) - z(0)., 
x [tanh( K \l - { 7 )+tanh( - v '- )], (11) 



2k B T 



2k B T 



where z(9) = [j,B+(vFq/2) cos(0— q ), and 6> q is the polar 
angle of the pairing momentum q. In order to determine 
the (second order) phase boundary between the normal 
and FFLO states, one needs to find the solution of ( |Tl| ) 
with A = and the largest possible B. At T = 0, we 
find the initial pairing instability occurs at q — 2/iB/vf, 
9q = 0,7r or ±7r/2 (i.e., for q along directions of maxi- 
mum gap), and [iB ss 1.06 Aoo- At lower fields the sys- 
tem undergoes a first order transition to the q — state. 
As the condensation energy of the FFLO state is quite 
small (the state is gapless over much of the Fermi sur- 
face), an excellent estimate for this field is simply the 
value obtained earlier for the level crossing between the 
normal state and the q = state, which leads to the 
window 0.56Aoo ^ iiB < I.O6A00 for the stability of the 
FFLO state, which is considerably larger than the win- 
dow A/V2 < fiB < A for the 2D s-wave case @. At 
finite temperatures, we find for T/T c < 0.06, the direc- 
tion of pairing momentum q remains the same as that of 
T = 0; however at T/T c > 0.06 this changes discontinu- 
ously to q = ±7r/4 or ±37r/4, i.e., q now points in the 
directions of minimum gap. At finite T we again use the 
crossing point of the free energies of the normal and q = 
paring states to estimate the boundary between the zero 
and finite momentum pairing states (Fig. ^), this window 
narrows and the magnitude of the pairing wavevector for 
the high field instability decreases and approaches zero 
at T = 0.56T C , where the high field phase boundary and 
the coexistence line between the q — state and the nor- 
mal state come together. In the FFLO phase, there is 
presumably a first order phase boundary across which 
the direction of the pairing momentum changes, which 
ends at the phase boundary separating the normal and 
FFLO states at T/T c w 0.06. In the present work we 
have not attempted to study this phase boundary. For 
T > 0.56T" C , there is no region with q ^ pairing and 
there is a continuous transition directly from the normal 
state to the q — state. This topology is also identical 
with that in the s-wave problem. We note that our re- 
sults on the phase boundary separating the normal state 
and FFLO state agree with a previous study by Maki 
and Won |15|, which are also confirmed in more recent 
work (Ref . |l6f| ) ; however in these works no estimate was 
given for the phase boundary separating the FFLO state 
and the usual zero momentum pairing BCS state. 
Application to the cuprates: Our analysis has been 
purely two dimensional, and for such systems a magnetic 
field parallel to the plane would behave precisely as a Zee- 
man field. For layered systems the situation, analyzed in 
some detail by Klemm et al. JT^] , is more complicated 
for orbital effects become important near T c where the 
interplanar coherence length is large. However, in the 
same approximation, there exists a lower temperature 
T* , where the vortex cores fit between planes and the 
orbital H C 2 diverges. Below T* (estimated as 0.84T C and 
0.99T C for YBCO and BSCCO respectively pi), assum- 
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ing two dimensionality should be an excellent approx- 
imation and hence the Zeeman response should domi- 
nate. Another limitation of our analysis is the neglect 
of scalar impurity scattering which destroys the FFLO 
state in dirty superconductors, and of spin-orbit scatter- 
ing which attenuates the pair-breaking effect of the Zee- 
man field. On these fronts, the news seems to be good: 
the state of the art YBCO and BSCCO samples are in 
the clean limit, and their residual resistivities translate 
into scattering times of order r ~ 10~ 12 s, which lead via 
the Elliott estimate (T3|, to a spin-orbit scattering time 
T SO ~ r/(Ag) 2 ~ 10"~l°s (Ag = g - 2 « 0.1). Conse- 
quently, tso > h/ {q^bB) for fields above a tesla and the 
neglect of spin-orbit scattering should not be too serious. 
A final caveat is the conventional BCS weak-coupling na- 
ture of our analysis, which is evidently problematic in the 
cuprates; absent a solution of the larger problems in the 
field, we are unable to say very much more on the issue. 

Nevertheless, the qualitative physics uncovered by our 
model calculation, should be quite robust to any mecha- 
nism that yields a d x 2_ y 2 state in a single layer |pof . Ex- 
perimentally, the low field effects discussed here should 
be readily observable, e.g. we estimate an enhanced spe- 
cific heat of magnitude Q.IHT mJ/mol-if 2 (k B T < fiB) 
from the existing data on YBCO 21 1, while the high field 
phase transitions and the FFLO phase would appear only 
at fields of order 100T, which are currently out of reach. 
Finally, while we have concentrated entirely on the paral- 
lel geometry, it is clear that a full account of the response 
at arbitrary orientations of the magnetic field will need 
to take account of the Zeeman physics discussed here. 
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FIG. 1. Fermi surface pockets (shaded regions) produced 
by a Zeeman magnetic field. The electrons in the pockets are 
unpaired and spin polarized along the direction of the field. 
The dashed lines indicate the extent of the smearing of the 
Fermi surface by the superconducting order at zero field, and 
show that the lateral extrema of the "normal" pockets are 
bracketed by regions of paired electrons. 



FIG. 2. The temperature-magnetic field phase diagram for 
a d x 2_ y 2 superconductor. The solid line is the second order 
phase boundary separating the normal state and the super- 
conducting state. Above T/T c — 0.56, the superconducting 
state is the zero momentum pairing state while below it is a 
finite momentum pairing state at high fields, which gives way 
to the zero momentum pairing state by a first order transi- 
tion along the dashed phase boundary. At T/T c » 0.06, the 
direction of the pairing momentum at the phase boundary 
changes discontinuously from that of the gap maxima to gap 
minima; the kink in the phase boundary reflects this change. 
The lower (upper) dotted lines are metastability lines above 
(below) which the normal (zero momentum pairing) states 
become local minima of the free energy. 
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Figure 2 




